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Abstract 

It is suggested a topological hierarchical classification of the infinite many Localized 
phases figuring in the phase diagram of the Harper equation for anisotropy parameter 
e versus Energy E with irrational magnetic flux uj. It is also proposed a rule that 
explain the fractal structure of the phase diagram. Among many other applications, 
this system is equivalent to the Semi-classical problem of Bloch electrons in a uniform 
magnetic field, the Azbel-Hofstadter model, where the discrete magnetic translations 
operators constitute the quantum algebra £/q(sZ 2 ) with q 2 = e l2nuJ . The magnetic 
flux is taken to be the golden mean uj* = (\/5 — l)/2 and is obtained by successive 
rational approximants uj m = F m _i/F m with F m given by the Fibonacci sequence 
F m .[OUTP-00-08S, cond-mat/0011396] 



Keywords: Topological Classification, Harper, Azbel-Hofstadter, Phase Diagram, Fractal, 
Quantum Groups 

PACS: 03.65.Sq, 05.45.+b, 71.23.Ans 



1 Adrss. for original paper: Dept. of Physics, Th. Phys., University of Oxford, Oxford 0X1 3NP, U.K. 
Curr. Adrss.: Dept. de Fis., Univ. da Beira Interior, Av. Marques D'Avila e Bolama, 6201-001 Covilha 
Curr. Adrss.: CENTRA, 1ST, Av. Rovisco Pais, 1049-149 Lisboa, Portugal 



1 Introduction 



In this article we study the phase diagram E x e of the Harper equation [2] 

ip n+1 + V> n -i + 2e cos(27T0 n )^„ = Eip n , 

(1) 

0n = 4>o + ujn . 

The Harper equation is a Schrodinger operator of two discrete variables ip and and, among 
other application, is associated with the quantum algebra of discrete magnetic translations 
[/q(s/ 2 ). Namely it is obtained in the Azbel-Hofstadter model [3J E] of Bloch electrons 
(electrons on a 2-dimensional periodic potential) in the presence of a uniform magnetic field 
in the Landau gauge. Here the anisotropy parameter e is the ratio between the spacing 
in the two directions of the underlying 2-dimensional periodic potential lattice, u is the 
magnetic flux per unit cell of the potential lattice and the quantum deformation parameter 

i S g2 = gi2™ j3J[7|_ 

The Harper equation is second order on ip. In order to simplify the numerical treatment of 
this equation we can cast it into a first order equation by considering the transformation 
of variables x n = ip n -i/ip n - The resulting equation is equivalent to the Harper equation (JTJ) 
and is known as the Harper map [S] 
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x n - E + 2e cos(27T0 n ) ' (2) 
<fin = 0o + un . 

In order to study the phase diagram E x e we need to identify the Localized and Extended 
phases. In order to do so we use Aubry and Andre study 0, they proved that the Lyapunov 
exponent A is proportional to the localization length, specifically A = —27. For the Harper 
map A < 0, such that for A = the phase is Extended and for A < the phase is Localized. 
The Lyapunov exponent corresponding to the <fi dynamics is always and the one concerning 
to the x dynamics is given by 

1 N 

A = (3) 

n=0 

where y n = log = hgx 2 n+v 

There exists already in the literature classifications, both for the Extended phases and 
for the Localized phases ^l] of equations (0) and In the remaining of this article 
we propose a new hierarchical topological classification of the several Localized phases for 
irrational flux to* = (y5 — l)/2 based on the topological winding number of the attractors 
on the space of variables (0, x) . In order to study the irrational phase diagram we take 
successive rational approximants u m . The fractal structure of the phase diagram emerge 
naturally in this framework and we suggest a rule that reproduces it. As we will show in 
detail, the great advantage of this new classification is that its labels correspond to the 
magnitude of the several different phases. Here the phase magnitude is understood as the 
area of each different phase on the full phase diagram. Also the area (hence magnitude) of 
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the several phases coincide with its order of appearance in the phase diagrams for successive 
approximants uj m to uj* . 

This article is organized as follows, in section 2 we introduce the irrational phase diagram 
and successive approximants as well as the one-dimensional attractor in the cf) x x diagram 
for irrational uj* . In section 3 we justify the identification of the <fi x x with a torus and 
introduce our classification. In section 4 we analyse the phase diagram and give the rules 
to build the Hierarchical Tree for the phase diagram. In order to better explain the fractal 
structure and how to build this Tree we have three appendixes. In appendix El we present 
the Lyapunov exponents for uj = uj* and e = 1 with enough resolution to identify the several 
phases (that correspond to each negative bump). In appendix El we list the x x diagrams 
up to 7] = 55 (with indication of the respective energies). Finally in appendix O we show 
graphically how to build the triangular structure up to level 9. 



2 Phase Diagram 

We proceed to present the numerical results and study the properties of the E x e phase 
diagram of the Harper map (J2J) which is equivalent to the Harper equation We consider 
successive rational approximants uj = uj m of the irrational value uj = uj* given by the golden 
mean. The rational succession uj m = F m _i/F m converges to the irrational golden mean 
uj* = (y/E — l)/2. Being F = F 1 = 1 and F m = F m _i + F m _ 2 the m th Fibonacci number in 
the Fibonacci sequence. 

Computing the phase diagram for successive rational uj m allows us to visualize the fractal 




Figure 1: Phase diagram for rational uj for <j)Q = 0. The Localized phases (L) are presented, (a) W5 = 5/8; 
(b) w 6 = 8/13. 

and denser. In the limit of irrational u*, the lines cross at the critical value e = 1 and 
the Lyapunov exponents go to zero at these points of intersection [10J. Although for the 
rational values uj m the phase diagrams depend on $0, in the irrational limit uj* the phase 
diagram is not sensitive to 0o- Note also that for e > 1 the Lyapunov exponents become 
always negative, this means that the phase transitions are smooth without undergoing an 
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Figure 2: Fractalization of the phase diagram. Lyapunov exponents A for e = 1 and (j>o — 0. (a) 
uj 7 = 13/21; (b) lu s = 21/34; (c) lu 9 = 34/55; (d) wio = 55/89; (e) tu = lo* = (y/E - l)/2; (f) lu = to*, zoom 
of(e). 

Extended phase. See [TT] for further details in these discussion. As a final remark we note 
that the phase diagram is symmetric with respect to both axes E = and e = 0. Also 
we stress that from the mathematical point of view the parameters E and e can be both 
negative or positive, from a physical point of view they should be positive, E is the energy 
and e is the ratio of the unit cell sides. The phase diagram becomes in the limit of irrational 
uj* a fractal (see figure 121 and |3J) , more exactly, a Cantor set j^J. It is this particular limit 
that interest us. From the above discussion it is enough to restrict our analysis of the 
several phases to the first quadrant for e < 1 and 0o — 0. 

It remains to study the attractors in the space of variables (<ft, x). Concerning the purpose 
of this work we describe the behaviour of the map in the Localized phases only, for further 
details on extended phases see [HI EI- For a given oj m , the map has 0-dimensional attractors 
with period F m , this means that the attractor is a _F m -cycle. In other words the map 
iteration cycle trough F m distinct points. In the limit u — > uj* it becomes an one-dimensional 
attractor F^, it is a line. These results are present in figure |U It is this property that we 
are exploring next in order to define our phase classification. 



3 Topological Classification 

It is time to define the classification. We are considering a map from the generic space of 
variables x x to a torus such that the identifications : = 1 and x : +oo = — oo hold. 
We note that we don't actually have to impose these identifications, they naturally come 
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Figure 3: Phase diagram for ui*. The first 12 zones are labeled in the diagram, 
from the definitions of the variables and x. 

The identification is justified by noting that 2-7T0 is an angular phase, then the variable 
is compact in the interval [0, 1] with identification : = 1 as any usual angular variable. 

The x identification x : +oo = — oo is somehow more subtle. By inspecting the attractors 
we realize that the one-dimensional attractors cross from ±oo to =Foo at the same value of 
(as an example see figures 0] and EJ) . The x x diagram gives us local properties of the 
original continuous variable ifi at each spatial index n. We note however that we completely 
loose all the information about the spatial localization itself (ordering on the spatial indexes 
n) of where (and in which order) the original variable ip have those properties. We are going 
to show that the x n = 0V-i/Vv variable is describing exactly the same situation for both 
cases x = +oo and x = — oo. The setup here is the original continuous variable being null 
■0 = at some (or many) points, x n = ip n -i/ip n = ±oo means that ip n = 0. To show the 
validity of this statement based on geometrical arguments see figure |H1 where two points, at 
different spatial localizations n and n', are considered, the discretized points ipn-i, i> n , Vw-i 
and are represented in the figure. Take in both cases |^n-i| ^ \^n\ an d iVV-il 3> W n \- 
This means that x n ~ +oo and x' n ~ — oo do correspond to the same situation of the 
continuous variable ip = 0. Both kind of situations will always occur for irrational u (many 
times for large N) because the period of the variable ip and the period of the underlying 
discretized lattice (given by the magnetic flux to) are incommensurate pH 1^1 HUj. 

In order to formalize our argument we can invert the map x n = ip n -i/tp n between variables 
obtaining ip n = ip n _i/x n , then we have that x n — > ±oo correspond to one only value 
ip n — > 0. Although not in the scope of the present study, it is interesting to note that 
x n — > correspond to ip n — ► ±oo, therefore a divergence of the original variable ijj. 
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Figure 4: Attractors in the space of variables </> x x corresponding to the Localized phase 1, for 0o = 0, 
E = 0.75 and e = 0.5. (a) uj 5 = 5/8; (b) lu 6 = 8/13; (c) oj n = 89/144; (d) u = lo* = (VE - l)/2, in this 
limit the Fm -cycle becomes a line. 
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Figure 5: Maps for lo* and e = 1. (a) r) = 11, E = 2.342; (b)r] = 10, E = 2.347. 



So we just proved that the variable identifications : = 1 and x : +oo = — oo are 
natural with out need for any a priori assumptions. Then the space of variables <fi x x is 
topologically a torus. Therefore the one-dimensional attractors, for irrational u*, are close 
contours living on the surface of that torus. By inspection we conclude that the attractors 
wind several times along the x coordinate. Furthermore the winding number along the x 
homology cycle is unique for each continuous Localized phase. So we classify these several 
phases accordingly, such that our classification label 7] is the winding number along the x 
homology cycle of the one-dimensional attractors. Note that the winding number in the <fi 
homology cycle is always 1 except for E = and e < 1 where there are two one-dimensional 
attractors with winding number in the x axis and 1 in the <fi axis ^U]. This particular 
case E = is excluded from our classification. 

The great advantage of this classification is that it also provide us a natural hierarchy based 
on the areas of the several phases. It coincides with the ordering of the areas occupied by 
the corresponding Localized phases. Furthermore as we consider successive u; m 's, the order 
of appearance of the phases also coincides with their hierarchical importance. As explained 
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Figure 6: The wave function is null, ip = 0, at two distinct spatial points but the corresponding x have 
opposite signs. 

in section El this is due to the Localized phases emerging from the splitting of the Extended 
phases. The diagram for irrational u* and the classification of areas up to 77 = 12 is 
presented in figure El As an example for higher 77 attractors we present in figure El the cases 
for 77 = 10 and 77 = 11. 

4 Hierarchical Fractal Structure 

Here we will define the generic rules to build the fractal structure of the phase diagram in 
the first quadrant (E > and e > 0). 

First let us state some conclusions obtained by inspection (CI) of successive phase diagrams 
for approximants u) m (see figures [T] and EJ) : 

CI.l In each new diagram for successive approximants the last phase on the far right (for 
E > 2) has always rj = and correspond to the phase with greater area. It is 
always present and doesn't contribute anything to explain the structure. So r\ = is 
excluded from the classification. 

CI. 2 The first phase to appear next (corresponding to the next greater area) as we move 
in the Fibonacci succession, is 77 = 1. Although we have two Fibonacci elements 
corresponding to this number, we have only one phase in the phase diagram. In some 
way (that we will explain later) this means that 77 = 1 is a degenerate phase. 

CI. 3 We recall from section |21 that for each approximant u m there are E m _i Localized 
phases, from the Fibonacci succession we have that F m = F m _ 2 + F m _i, then as we 
go trough the approximants, from u m to u m+ i emerge exactly F m ^ 2 new Localized 
phases. It seems then a good guess to organize the phases in levels such that in each 
level k = m — 1 there are exactly F^_ 2 labels in the range F^-i < rjievei k < Fk- We 
are therefore building a triangular structure. 
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CI.4 Also as we increase m the last phase on the far right (just before 77 = 0) corresponds 
always to one label belonging to the Fibonacci succession 77^ = F m with k — m — 1. 
We will call r]J} the right pivots. 

CI. 5 The last phase on the left side of the diagram (near E = 0) only changes for the 
second odd number in the Fibonacci succession (even, odd, odd, even, odd, odd,. . .). 
This means that only for k = 3n + 1 (for all integer n) we will have a label f]^ n+1 = 
F 3n + F 3n _ 2 . We will call r)% the left pivots. 

CI. 6 As we move along the approximants new phases will emerge in the bulk of the diagram 
next to the ones that already exist in the previous approximants. However if one phase 
rj n appears at some approximant u m , at the immediately next approximant u m+ i no 
new phases will appear next to this phase rj n . So, for each new phase appearing, 
there are not adjacent new phases at the immediately next approximant. 

CI. 7 Also these new phases appearing adjacent to the already existing phases in the bulk 
will appear either on the right or on the left of the already phases that already exist. 
So the order of appearance of new phases alternate between left and right. 

So we have a triangular Hierarchical Tree organized in levels, each level k has exactelly Fk-2 
elements (see point CI. 3 above) and we are ready to define the rules to build it. In order to 
accomplish it we first give the rules for the edges of the triangle and then define a recursive 
rule for the bulk of the triangle. For the reader interested in the full details we present 
in appendix [X] the Lyapunov exponents for u = u* and e = 1 with enough resolution to 
identify the several phases (that correspond to each negative bump) and in appendix iBl we 
list the (f) x x diagrams up to rj = 55 with indication of the respective energies. 

We can define the vertex and edges of the triangle following the following rules: 

R.I On the top vertex of the triangle is placed the degenerate label rj = 1 (see point CI. 2 
in our list of conclusions above). It is degenerate simply because it belongs simulta- 
neously to the right and left edge of the triangle, i.e. it is both a left and a right pivot 
Vo = Vo = 1 ( see conclusions by inspection CI. 2, CI.4 and CI. 5). 

R.II At the right edge of the triangle we have at each level k the right pivots = F^+i 
distributed, in sequence (see point CI.4 on the list of conclusions above). 

R.III On the left edge not all levels have a pivot, we will have left pivots only for level (3n+l ) 
(for any integer n) with labels r] 3n+1 = F 3n + F 3n _ 2 (see figure EJ). 

Next we will define a recursive rule for the bulk of Tree. In order to distribute the labels 
in each level, consider the Tree completely built up to level k and let us construct the 
next level k+1. We have therefore Fk-i labels in the range < T]i eve i k+i < -ffc+i (see 
point CI. 3). Although we already have rules for the edges of the triangle (see Rule II and 
Rule III) we include in the recursive rules also the edges. The following rules should be 
applied in sequence: 

R.l Group the two last right pivots (Vk-uVk) an d sum them together, the result is the 
new pivot of level k + 1, i.e. rjj} +1 = rjj}^, r/j}. 
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R.2 Group in pairs all the remaining elements (the labels) up to level k — 1 (we note that 
no other elements of level k will be used except for the right pivot r]j}). The way to 
group them is to start from the right pivot T]^_ 2 and group it with the immediately 
adjacent element to the left. Then we take the next element immediately adjacent to 
the left and group it to next one adjacent to the left, and so on. Here adjacent means 
the order in the phase diagram as we move to the left, in terms of the Hierarchical 
Tree, i.e. the closest element in relation to the horizontal projected distance. From 
two in two levels there will be out of the pairs a single left pivot (see point CI. 5 above), 
then after applying Rule R.3 below, apply also Rule R.4. 

R.3 Sum the right pivot of level k (rjj}) to all the pairs obtained in Rule 2, then we swap the 
order of the elements inside each pair and place the resulting (swapped elements) pair 
at level k + 1 between the elements of the original pair. Again, here, between means 
in relation to the horizontal projection of the elements. We note that these rules do 
reproduce the behaviour described in the conclusions by inspection CI. 6 and CI. 7. 

R.4 From two in two levels we will have a single left pivot unpaired of level k-1, in this 
case we sum it to the right pivot rj^ and place the result as the left pivot of level k+1. 

We note that in order to apply the recursive rules we only need the two first labels, i.e. 
Vo = Vo = 1 an d Vi — 2- We note that Rule II is equivalent to Rule 1 and Rule III 
to Rule 4, we simple list both set for completeness. In this way the rules defining the 
triangular structure are complete! The first 9 levels of the Hierarchical Tree are shown in 
figure [3 we can compare it with the corresponding phases in the phase diagram (figure EJ) • 
In appendix Owe show graphically how to build the triangular structure up to level 9 using 
these recursive rules. 

5 Conclusions 

In this paper we have proposed an integer topological hierarchical classification for the 
infinite many Localized phases in the phase diagram Exeof the Harper map for irrational 
parameter us* = (V5 — l)/2. This map is equivalent to the Harper equation which, among 
other applications, describes the Azbel-Hofstader model with E being the energy, e the 
anisotropy parameter and u the magnetic flux per unit cell. The classification is based in 
the winding number around the holonomy cycles of the One-dimensional attractors in the 
toroidal space of variables <fi x x. For the classifications obtained, the fractal structure of 
the E x e phase diagram is built and organized into a triangular structure that exactly 
describes the phase diagram taking in account the ordering of the several phases. Moreover 
the integer classification labels (1, 2, 3, ... ) exactly correspond to the phase importance, 
concerning their area in the phase diagram and order of appearance in the succession of 
phase diagrams for the several approximants uo n converging to u*. So in this since the 
triangular fractal structure correctly describes all approximants. As possible future works 
in this subject it would be interesting to investigate the fractal structure for u given by 
other irrational numbers and, if possible, to find a rule for generic irrational numbers and 
respective approximants. 
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Figure 7: Fractal structure of the phase diagram (to be seen in landscape orientation). 
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A Lyapunov Exponent 



In this appendix we present the Lyapunov exponents for e = 1 and on*. The several phases 
are labelled until r\ = 55. 
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Figure 8: Lyapunov exponent A as a function of energy E for e = 1 and = a;*. The 
Details I to V correspond to the marked regions in the first graph. 



B (j) x x diagrams 

Next we list all the attractors in the space of variables x x up to rj = 55. 
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Figure 9: Attractors in the space of variables x x up to 77 = 21. 
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C Building the Hierarchical Tree 



We show graphically how to build the Hierarchical Tree up to level 9. Level 1 and Level 2 
are given and only contain the vertex (simultaneously left and right) pivot t)q = vfc = F\ = 1 
and the right pivot r]f = F 2 = 2. We give details of the application of the construction 
rules up to level 6 which contains all rules. The remaining levels do not involve any new 
rule and can be build straight forwardly. 

On the following figures the pivot used to construct each level n + 1 (i.e. 77^, the pivot 
of the level immediately above the one we are constructing) is inside a darker circle. The 
remaining elements used are grouped in pairs pictorially marked. We note that up to some 
level n we have present in our Tree all labels up to rj^ = F n+1 . 



level 1 



level 2 



level 3 




Figure 12: Building level 2 (n+ 1 = 2). The right pivot at level n = 1 is r?f = F 2 = 2. The 
next pivot is obtained by adding the two previous pivots = r/^+r/f = F ± +F 2 = 1 + 2 = 3. 
There are no more elements in this level. 
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level 2 



level 3 



level 4 




Figure 13: Building level 3 (n + 1 = 3). The pivot at level n = 2 is = F 3 = 3. The next 
pivot is obtained by adding the two previous pivots rj^ — rjl + r\ 2 — 2 + 3 = 5. There are no 
enough elements to group in pairs, we have one single element at the vertex of the triangle 
t?q = t)q = Fi = 1 that we deal as being a left pivot. Then we sum it the right pivot of 
level 2 r\2 obtaining the new left pivot 773 = 770 + — F x + F 3 — 1 + 3 = 4. 
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level I 
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Figure 14: Building level 4 (n + 1 = 4). The pivot at level n — 3 is r]^ — F4 — 5. The 
next pivot is obtained by adding the two previous pivots rjf — T]\ + 773 = 3 + 5 = 8. The 
remaining elements on the levels above (not including) level 3 are grouped in pairs, we 
obtain one single pair (1, 2). We sum this pair to the current level pivot rj§ = 5 and place 
the resulting pair with the order of the elements reversed below and between the elements 
of the original pair (1,2). This means we obtain the pair (7, 6) and place it as shown in the 
figure. 
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Figure 15: Building level 5 (n + 1 = 5). The pivot at level n — 4 is 77^ = F 5 — 8. The 
next pivot is obtained by adding the two previous pivots 77^ = 773 + r\\ = 5 + 8 = 13. The 
remaining elements on the levels above (not including) level 4 are grouped in pairs, we 
obtain two pairs (4,1) and (2,3). We sum these pairs to the current level pivot 77^ = 8 
and place the resulting pairs with the order of the elements reversed below and between 
the elements of the original pairs. This means we obtain the pairs (9, 12) and (11, 10) and 
place them as shown in the figure. 
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Figure 16: Building level 6 (n + 1 = 6). The pivot at level n — 5 is r)^ = F 6 — 13. The 
next pivot is obtained by adding the two previous pivots rj^ = i]^ + r}\ = 8 + 13 = 21. The 
remaining elements on the levels above (not including) level 5 are grouped in pairs from 
right to left. As expected for level 6 we have one lonely left pivot i]^ = 4 and we obtain 
three pairs (1,7), (6, 2) and (3, 5). We sum the current level pivot r/^ = 13 to the lonely left 
pivot obtaining the new left pivot = V5 + V3 = 13 + 4 = 17. As for the three pairs we 
sum them to the current level pivot 77^ = 13 and place the resulting pairs with the order 
of the elements reversed below and between the elements of the original pairs. This means 
we obtain the pairs (20, 14), (15, 19) and (18, 16) and place them as shown in the figure (to 
be seen in landscape orientation). 
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Figure 17: Building level 7 (to be seen in landscape orientation). 
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Figure 18: Building level 8 (to be seen in landscape orientation). 
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